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We investigate recently proposed method for locating critical temperatures and introduce some modifications 
which allow to formulate exact criterion for any self-dual model. We apply the modified method for the Ashkin- 
-Teller model and show that the exact result for a critical temperature is reproduced. We test also a two-layer 
Ising model for the presence of eventual self-duality. 



1. Introduction 

The recent proposition of the analytical 
method for locating critical temperatures in some 
spin systems employs moments of the transfer 
matrix Let us remind briefly its basic as- 
sumptions. 

For a d-dimensional spin system we define the 
characteristic function of rank n 



(1) 



where T is a transfer matrix of our system with 
linear size L. The latter quantity was introduced 
for regularization purposes. The dependence on 
the inverse temperature (3 = 1 /T is hidden in the 
definition of the transfer matrix T . The same 
method relies on the hypothesis that the location 
of a maximum of the function p occurs at a critical 
point ( "maximum criterion" ) 

Pc Pmax- (2) 

It remains to test this very attractive hypothe- 
sis in practice. Among others an important ques- 
tion is: how big is the class of models for which 
the relation (||) indeed holds for any n ?. A con- 
jecture was given, that a maximum rule is at least 
valid for self-dual systems. Indeed it was checked 
directly for simple self-dual systems like Ising and 
q state Potts models on a square lattice ||. 

Recently Souza et al. have analyzed the 
Ashkin- Teller model and found that the criterion 



does not reproduce correctly its phase structure 
||]. However in particular case of the isotropic 
Ashkin- Teller, where the phase structure is de- 
termined from the self-duality condition, the pre- 
dicted critical temperature is a good numerical 
estimation of the exact value. 

In this paper we introduce some modifications 
which allow to formulate exact criterion for any 
self-dual model. We show it on the Ashkin- Teller 
model as an example. Apart from determining 
critical temperatures our method can be used to 
test the existence of duality transformations for 
more complicated models. We investigate a cer- 
tain kind of duality relation in the two-layer Ising 
model and show that it is not self-dual. 

2. Application to the Ashkin- Teller model 



The Ashkin- Teller model 
Hamiltonian 



is defined by the 



H = -J - J / J O-jCTj - J' / J SiSj 

(id) 

(id) 



(3) 
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with two Ising-like variables Si and <r,; located on 
the same site i of a quadratic lattice. In some 
sense the AT model interpolates between the Ising 
(J4 = 0) and 4-state Potts model (J = J' = J4). 
In general its phase structure is rather compli- 
cated and can not be simply obtained. 

In the simpler isotropic AT model ( J = J') the 
duality transformation provides for J > J4 the 
position of a critical point. It is convenient to 
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introduce dimensionless coefficients 

K = /3J , K = /3J, K4 = f3J 4 

and the quantities 

ujq = exp[2K + K A + K ), 

lu 2 = exp (—Ki + K ), 

lj 3 = exp(-2K + K 4 + K ). 



(4) 



(5) 



Then the partition function for the isotropic AT 
model is invariant 



under the duality transformation 
uj 3 = ~(u>o - 2uj 2 + ai 3 ). 



(6) 



(7) 



The critical temperature is known from the self- 
duality condition 



uiq = 2uj 2 + ^3- 

From definition (Q) follows that 
i 

Pn (f3) = lmi ( |£ 



(8) 



(9) 



where Z\ and Z n are the partition functions of a 
chain of length L and of such n coupled chains 
respectively. In our example the characteristic 
function p depends on two couplings K and K4, 
(as can be easily noticed Kq does not appear in 
the final result). Therefore we may expect similar 
to ffl) invariance for the function p 



p n (K',K A )=p n (K,K 4 ) 



(10) 



where K' and K 4 are the dual couplings to K and 
K4 respectively. 

For simpler models (like Ising or Potts) the 
function p depends only on one coupling or equiv- 
alently on (3. Then it can be easily shown that if 
there is only one maximum of p at some (3 m ax it 
must coincide with the self-dual point (3 max = (3*. 
The same is not true in our case. In fact the 
straight line parametrized by 



is not consistent with the duality transforma- 
tion i.e. the dual image of a point (K, K4) taken 
from line ([□]) does not lie necessarily on this 
line. Our proposition of the modified criterion 
employs the duality transformation to find the 
proper direction for passing the self-dual curve. 
Since we are mainly interested in the behav- 
ior of the function p near the real critical point 
it is sufficient to consider the linearized dual- 
ity transformation in a vicinity of some self-dual 
point (Kq, K* 7 KX). Let us denote the vector 
k T = (K -K£,K- K*, if 4 - KX) and its image 
under the duality transformation as k' . Then 



k' = Mk 

where the transformation matrix M is: 



-a 2 







'b-b 




(12) 



(13) 



with a = sinh2iir* and b — cosh2if*. The matrix 
M has three eigenvectors: two with eigenvalue 1 
and one with —I. The latter eigenvector 



/ 2cosh2.Fi:* sinhTftT* 



V cosh TFT* - 3' 3 - cosh TFT* 



1 



(14) 



is interesting since it corresponds to approaching 
the self-dual curve from both sides and provides 
the valid parametrization (analogous to (|lT|) ) in 
the modified criterion 
3-cosh4X* 



K 4 
K 



sinh AK* 



-(/3-K*)+KZ, 



0, 



(15) 



K = (3 J, Ki = f3J 4 



(11) 



where J = I is assumed. 

Of course the arguments given above are not 
a rigorous mathematical proof. In particular the 
assumption that for modified criterion the func- 
tion p has only one maximum has to be checked 
directly for a given model. For the AT model few 
lowest moments can be easily computed numeri- 
cally. Table. I contains the critical temperatures 
obtained from the p2 function in the original and 
modified method. In accord to Q original cri- 
terion is exact only for J4 = (two decoupled 
Ising models) and for J4 = I (four state Potts 
model). The modified method gives exact results 
for J4 < I up to rounding errors of order of I0~ 8 . 
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J 4 


ftcrit 


ftrnax 




0.00 


0.44068679 


0.44068679 


0.44068679 


0.10 


0.41215166 


0.41150224 


0.41215165 


0.20 


0.38799451 


0.38716551 


0.38799451 


0.30 


0.36721416 


0.36642983 


0.36721416 


0.40 


0.34910080 


0.34845723 


0.34910080 


0.50 


0.33313596 


0.33266131 


0.33313596 


0.60 


0.31893130 


0.31861775 


0.31893130 


0.70 


0.30618959 


0.30601070 


0.30618958 


0.80 


0.29467908 


0.29459934 


0.29467908 


0.90 


0.28421612 


0.28419628 


0.28421611 


1.00 


0.27465307 


0.27465307 


0.27465307 



Table 1 

Comparison of different estimations of the crit- 
ical temperature for the isotropic Ashkin- Teller 
with J = 1: (3 cr it (exact value) (3 max (original 
criterion), $ max (modified criterion). 



For J4 > 1 the self-dual curve (g) is not longer 
a critical line. In fact the phase structure is more 
complicated with two second order phase tran- 
sitions. Our criterion predicts the existence of 
only one phase transition exactly on the self-dual 
curve. 

3. Application to the two-layer Ising model 

It is very interesting to apply the modified cri- 
terion to other models (not necessarily self-dual) 
for which the original maximum rule does not give 
exact answers. The two-layer Ising model is an 
important example. Its Hamiltonian resembles 
that for the AT model 



H = — J <7j<Tj 



J SiSj 



J 2 > Siffi 



(16) 



The main difference is in the form of the interac- 
tion term and in the interpretation of spin vari- 
ables s and a which are now assigned to two dis- 
tinct layers. Neither the duality transformation 
nor the exact critical temperature is known for 
this model. One can still attempt to determine 
the slope of the straight line in the modified cri- 
terion from the condition 



Per it 


n=2 


n=3 


n=4 


n=5 


0.2755 


0.6399 


0.5969 


0.5915 


0.5929 


0.2760 


0.6233 


0.5746 


0.5647 


0.5625 


0.2765 


0.6067 


0.5526 


0.5382 


0.5324 



Per it ~ P 'n 



(17) 



Table 2 

Slopes of the lines (see eq. (|i~5|)) in the modified 
criterion for the two-layer Ising model. 



A priori, it is not even certain wheather the re- 
lation ( |j~7| ) can be satisfied at all. We found 
that we were able to estimate the slope numer- 
ically for few lowest moments (Table. 2). For P cr u 
( J = J4 = 1) we take the value 0.2760, the result 
of Monte Carlo simulations 0]. Since it is not 
exact, with the statistical errors on the last deci- 
mal digit, for comparison the results for two other 
possible values of (3 cr u are also given. The result- 
ing slopes depend on n which is exactly what we 
expect for not self-dual models. Since we do not 
have at our disposal sufficiently precise estima- 
tions of the real critical temperature it is prema- 
ture to present similar values for other critical 
points [K, K^)- However it is rather obvious that 
the results would not differ qualitatively. 

In conclusion, the method of moments is very 
closely related to self-duality of a given model. 
For any self-dual model there is a unique direc- 
tion in the space of parameters which allows to 
construct exact criterion. Conversely if the direc- 
tion does not depend on n then this is a strong 
argument for the existence of a self-duality rela- 
tion in a model. 
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